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Abstract 

A method for measuring translational diffusion constant in liquids via NOON state 
is described using a quantum circuit and is experimentally demonstrated using a 
model system. When compared with the standard single quantum method, the 
NOON state method requires shorter diffusion delays and weaker gradients. These 
improvements depend on the on the size of the NOON state. Due to the high 
sensitivity of the NOON state for the changes in the local magnetic fields, this 
method enables studying slow diffusion and studying diffusion with limited strengths 
of pulsed-field-gradients. 
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1 Introduction 



Driven by the internal thermal energy, the atoms or molecues of a bulk matter 
exhibit random translational motion, which is termed as translational diffu- 
sion [TP] . Diffusion explains the net flux of particles moving from a region 
of higher concentration to the lower concentration, although diffusion process 
exists even when there is no macroscopic concentration gradient. Diffusion is 
a fundmental transport mechanism in liquids and gasses and therefore mea- 
surement of diffusion constant is important for understanding many physical, 
chemical, and biochemical processes [SBjlSllG] . The diffusion constant {D) is de- 
scribed as the amount of a particular substance that diffuses across a unit area 
in unit time under the influence of a unit concentration gradient [1] . According 
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Fig. 1. Hahn-echo pulse sequence with two equal PFGs for the measurement of 
diffusion constant. 

to the Stokes-Einstein theory, for a fluid with visocity rj aX a. temperature T, 
the diffusion constant is given by 



where k is the Boltzman constant and is the Stokes radius of the diffusing 
particle [TIH] . The denominator in the above expression is termed as the friction 
coefficient. 



Diffusion constant of a liquid can be measured by NMR either with the help 
of relaxation studies or more conveniently using pulsed-field-gradients (PFGs) 
[9] . Fig. [1] shows a standard pulse sequence based on the Hahn-echo method 
for measuring diffusion constant [TOlITT] . In this method, maximum amplitude 
of Hahn-echo can be obtained by a refocussing vr pulse placed between two 
PFGs if the molecues do not change their place during the period A. Dif- 
fusion renders the molecules to change their place and therefore for a fixed 
duration A, the echo amphtude decreases with the strength of the PFGs. Such 
PFG methods are widely used for the measurement of diffusion constants and 
already numerous improved sequences are available [T2|13fl4] . For example, 
replacing Hahn-echo by stimulated-echo sequence reduces the decay due to 
transverse relaxation and results in enhanced echo amplitudes [TSflG] . Appli- 
cation of long-lived singlet states allows one to study slow diffusion requiring 
long intervals (A) between the PFGs [T7]. Fast single-scan measurement of dif- 
fusion has also been demonstrated by realizing z-coordinate dependent PFG 
strengths using frequency swept refocusing pulses [T8] . 

In this letter we describe another approach inspired by quantum information 
theory to speed-up the measurement of diffusion constant. We propose to pre- 
pare the diffusing spin system into a NOON state. Although this method has 
similarities with the multiple-quantum method for studying diffusion p^20f21] 
the current method differs in actual implementation and its scope. We briefly 
outline the theory in section 3 and we describe the experimental demonstration 
in section 4. 
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2 Diffusion constant by NOON state 



2.1 NOON state 



We shall follow the notations of quantum information theory where ±1/2 
states of a spin 1/2 nuclues is represented by basis kets |0) and |1) of a quantum 
bit or a 'qubit'. The NOON state of an A^-qubit system is a superposition of 
all the N qubits being in state with all those in state 1: 



\NOON) = {\N, 0) + |0, N))/V2 
= (|00---0) + I11---1 



[22]. Recently several applications of NOON states have been discovered [23|24f25] . 
It has also been used for sensing weak magnetic fields |26j . 

The circuit for the preparation of NOON state, shown in the first part of Fig. 
[21 consists of a Hadamard gate (H) and a CNOT gate. This circuit acts on 
a quantum register with a single 'control' qubit and a set of (A^ — 1) 'target' 
qubits initialized in |00 ■ ■ - 0) pure state 

|00 ■ ■ -0) ^(|0) + |1))|00 ■ ■ ■ 0) 

v2 

CNOT 1 Q00---0) + |ll---l)) = \NOON). (3) 
v2 



Preparing a pure NOON state is out of bounds from a highly mixed thermal 
equilibrium state that is available in liquid state NMR systems at ordinary 
temperatures. However, it is rather easy to prepare a pseudopure NOON state 
which is of the form 



PmoN = + e\NOON){NOON\, (4) 

where the scalar quantity e is a measure of purity, which is of the order of 10~^ 
in typical NMR setups [27]. The identity part does not evolve under unitary 
dynamics neither does it give raise to NMR signals. Therefore we ignore this 
part as a background and continue using ket notation |iV00iV) for simplicity. 

The CNOT implementation is particularly easier with a heteronuclear spin sys- 
tem of type AMjv-i, wherein each M spin of the A^ — 1 magnetically equivalent 
spins, is having an indirect spin-spin interaction with the common heteronu- 
clear spin A with a coupling constant of, say J Hz. We shall use the single 
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Fig. 2. Circuit for the measurement of diffusion constant using NOON state. NOON 
state is initialized using a Hadamard gate (H) and a CNOT (©) gate. Tlie diffusion of 
particles is monitored using two identical PFGs of strength Gi placed on either side 
of TT pulses. The second CNOT converts the evolved NOON state into an observable 
single quantum coherence of the 1st qubit. The signal of the first qubit is measured 
after the PFGs G2 and G3 select the coherence transfered from the NOON state. 

heteronucles as the 'control' qubit and the rest as the 'target'. Such a system 
allows parallel implementation of CNOT gate (upto a global phase) using the 
sequence 



t^CNOT = YiYuUjYl^UjYAXAYAYMXM (5) 

[28|29] . Here X = e"*^''/^)'^^ and X = e+*(''/^)'^^ are rotations on the spins 
indicated by their subscripts (/ is the spin operator and Y, Y are defined 
similarly). Under Uj = e"**^'^/^^'^^^ J2m^^ all the target qubits simultaneously 
undergo J-evolution with the control and is realized simply by a 1/ (4 J) delay. 
In our experiments we have replaced the Hadamard gate on A with a pseudo- 
Hadamard gate which can be implemented by a single Ya pulse. The rest of 
the circuit in Fig. |2] is described in the following. 



2.2 Phase encoding and dijfusion measurement 



Consider first the simple Hahn-echo sequence for measuring the diffusion con- 
stant as illustrated in Fig. [1] For simplicity, we ignore decoherence and other 
experimental errors in the following analysis. The first gradient encodes the 
longitudinal distribution of the sample such that a spin in position z gets a 
phase (f){z) = 'yzGzS, where 7 is the gyromagnetic ratio of the spin and 6 
is the effective duration (duration x shape-factor) of the gradient Gz- Un- 
der this gradient the pseudopure coherence |0) -|- |1) (ignoring the backgroud) 
evolves to |0) + e*''^|l). After the central vr pulse flips the basis states the second 
gradient encodes another identical phase so that the resulting state becomes 
6**^1 1) + 6**^10) = |0) + |1), up to a global phase. Thus if the spin does not 
change the position there is no change in the relative phase. However if, dur- 
ing the period A, the spin changes the position by an effective amount dz, the 
state after the second gradient will be |0) + e*"^"^!!), where the relative phase 
is d4> = ■jdzGzS. 
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Due to the random molecular motion of a large number of spins each of which 
acquiring a different relative phase, the average Hahn-echo does not acquire a 
net phase but diminishes in amplitude. Theoretically, the decay of the Hahn- 
echo signal S can be given as. 



S{G,) = Soexp{--f'GyD{A - 6/3)} , (6) 

where D is the diffusion constant, Sq = S{0) is the normalization factor, and 
A — (5/3 is the correction to A due to the finite durations 6 of the PFGs |12j . 

Similar analysis can now be carried out for the NOON state using the circuit 
shown in Fig. HJ Under the Gi — it — Gi sequence, a spin system in NOON 
state diffusing through a distance dz acquires a net relative phase and becomes 
(|00 ■ ■ ■ 0) + e^'^'^lll ■ ■ ■ The relative phase acquired is d(f) = -f^sdzGJ, 

where 7efr = {1 + (iV — 1)jm/ia}ia = l^A- The factor I is also known as 
'loopsidedness' of the NOON state. Larger the value of /, more sensitive is the 
NOON state for the phase encoding. Since NOON state is a multiple quantum 
coherence it is necessary to convert it back to single quantum coherence before 
detection. This conversion can efficiently be carried out using a second CNOT 
gate: 



^(|00 ■ ■ -0) + e'^^lll ■ ■ ■ 1)) ^(|0) + e''''^|l))|00 • • -0). (7) 
V 2 V 2 

Thus the phase encoding due to the diffusion, i.e., dcj) has been transfered to 
one transition of the control spin. However as explained before, for a large num- 
ber of molecules undergoing diffusion, the above phase encoding results in the 
attenuation of the control transition. Instead of starting with an initial pseu- 
dopure state, it is rather convenient to start with a thermal equilibrium state 
and use two PFGs G2 and G3 to select out the desired coherence pathway that 
is passing through the NOON state. The ratio of these two PFGs are adjusted 
depending on the relative gyromagnetic ratio: G3 = G2{{N — l)7Af/7A + !)• 

In the case of large loopsidedness /, the NOON state helps us to study diffusion 
with weaker PFGs and smaller durations (A) between them. An experimental 
demonstration of this method for a model system is described in the following 
section. 



3 Experiment 



The sample consisted of 100 [il of trimethylphosphite (see inset of Fig. [3]) 
dissolved in 500 ^1 of dimethyl sulphoxide-D6. All the experiments were carried 
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Fig. 3. The P spectra corresponding to single quantum excitation from thermal 
equilibrium (upper trace) and after converting the NOON state into the single quan- 
tum coherence using a CNOT gate (lower trace). The inset displays the molecular 
structure of trimethylphosphite (P(OCH3)3). 

out on a 500 MHz Bruker NMR spectrometer at an ambient temperature of 
300 K. Each of the nine magnetically equivalent spins are coupled to the '^^P 
spin via indirect spin-spin interaction with a coupling constant of J = 11 Hz 
thus forming an AMg spin system. Accordingly phosphorous spectrum splits 
into 10 lines as shown in Fig. |3l An initial INEPT transfer was used to enhance 
the ^^P polarization. Then, as described in the previous section, NOON state 
was prepared with a pseudo-Hadamard gate and a subsequent CNOT gate, 
and is converted back to a single quantum coherence using a second CNOT. 
The selection PFGs were adjusted such that G3/G2 = d'-yn/lp + 1 = 23.23, 
to select the 10-quantum coherence pathway. A two-step phase cycle of Xm 
pulse of the first CNOT along with the receiver phase helped to reduce artifact 
signals. The single transition in '^^P spectrum of Fig. [3] indicates the selection 
of NOON-state coherence pathway. 

We first carried out the single quantum diffusion experiment as described in 
Fig. [H The parameters of the experiment are shown in Fig. H^. While keeping 
all other parameters fixed, we varied only the strength of the PFGs Gi in 
the range to 0.3325 T/m in 20 equal intervals. The integrated intensity of 
the ^H signal as a function of Gi is shown in Fig. H^. The diffusion constant, 
D = (6.24 ± 0.06)10^^° m^ s~^ was obtained by fitting these data points to 
the Gaussian expression given in ([H]). 

Then we prepared the NOON state as explained above, and diffusion exper- 
iment was carried out as in Fig. O The main parameters of the experiment 
are shown in Fig. |4)d. The lopsidedness for this system is / = 9.4 and ac- 
cordingly the NOON state experiment requires much shorter A and 6 values 
to reach similar attenuation of the echo signal as that of the single quantum 
experiment. In FigHb, the integrated intensity of the single transition, after 
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Fig. 4. Intensity of the echo signals as a function of the gradient strength Gi with 
(a) single quantum coherence and (b) NOON state. The dots represent the 
experimental data and the lines represent the Gaussian fit. The insets in both the 
figures display the parameters of the experiments. 

the selection of NOON-state coherence pathway, is plotted against Gi. The 
Gaussian fit, again using the expression ([6]), with 7efr = I'Jh lead to the dif- 
fusion constant D = (6.17 ± 0.25)10~^° s~^. While the two methods gave 
identical values for the diffusion constant, the error-bar is slightly larger in 
the latter case due to the added complexity of the circuit, the imperfections 
in the selection of coherence pathway, reduced signal to noise, and due to the 
faster relaxation of the NOON state. On the otherhand, the shorter timescales 
of the NOON-state experiment suggests its possible applications in studying 
slow diffusion or studying diffusion with limited PFG strengths. 



4 Conclusions 

We described the application of NOON states in studying translational dif- 
fusion in liquids. Although the basic principle is general, the method is par- 
ticularly convenient in the presence of a set of magnetically equivalent nu- 
clei interacting with a heteronucleus via J-coupling. We have demonstrated 
the experimental measurement of diffusion constant in trimethylphosphite, 
which is a AMg spin system. Both single-quantum and the NOON-state ex- 
periments lead to identical values for the diffusion constant, but the error-bar 
was slightly larger in the latter case due to the additional complexities. How- 
ever, the NOON state experiment required weaker pair of PFGs and an order 
of magnitude shorter duration between them, indicating the possible applica- 
tions in studying slow diffusion. The other available method for measurement 
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of slow diffusion involves the use of long-lived singlet states [T7], which re- 
quires a pair of weakly interacting homonuclear spins spatially isolated from 
the rest of the spins [30j. In the current method we prepare a state that is 
highly sensitive to the gradient encoding of the z-coordinate and hence is able 
to capture the effects of diffusion even in shorter time-scales and with weaker 
PFGs. It might also be possible to combine the NOON state technique with 
the single scan techniques to achieve ultra-fast diffusion measurements. Such 
an experiment will then allow studying dynamic situations with time-varying 
diffusion constants. 
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